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Approximate  Analytic  Solutions  for  the  Primary 
Auroral  Electron  Flux  and  Related  Quantities 


I.  I NTRODl  CTION 

The  problem  of  calculating  the  distribution  function  for  auroral  electrons  may 
be  approached  from  two  points  of  view,  one  approach  is  to  neglect  collective 
effects  and  treat  short  range  particle-particle  collisions  using  linear  transport 
theory;  the  other  is  to  neglect  short  range  collisions  and  treat  collective  effects 
using  quasilinear  or  nonlinear  plasma  kinetic  theory.  At  high  altitudes  the  latter 
approach  is  often  used  as  the  plasma  is  essentially  collisionless,  whereas  at  low 
altitudes  {below  about  250  km)  linear  transport  theory  is  generally  used  as  the 
plasma  is  neutral  particle  collision  dominated.  The  question  of  how  to  combine 
the  short  range  and  collective  aspects  of  the  problem  is,  and  should  be,  an  active 
area  of  research.  In  this  paper  we  confine  our  attention  to  the  lower  ionosphere 
and  apply  the  methods  of  linear  transport  theory  to  calculate  the  auroral  primary 
electron  flux. 

There  is  an  extensive  literature  on  the  subject  of  elect ron  transport  in  the 
auroral  ionosphere.  The  various  approaches  may  be  categorized  as  somi- 
empirical,  range  theoretic,  Pokker-Planck,  Monte  Carlo,  and  transport  theoretic. 

Karly  work  on  transport  properties  of  KeV  auroral  electrons  comes  from 
Chamberlain.  Estimates  were  made  of  ionization  rates  as  a  function  of  altitude 

(Received  for  publication  2  March  1981 

1.  C  hamberlain,  J.  W.  (1961)  Physics  of  Aurora  and  Airglow,  Academic  Press,  N.  V 


for  energetic  KeV  electrons.  The  rates  were  based  on  deposition  results  of 

Spencer  who  solved  the  electron  transport  equation  using  the  continuous  energy 
9  3  4 

loss  approximation.  Rees  (a  semi-empirical  approach)  applied  an  energy 
dissipation  function  based  on  laboratory  data  by  Grun.  **  The  forms  of  the  func¬ 
tions  obtained  by  Rees  provided  altitude  profiles  of  the  energy  deposition  and 
ionization  rates  for  monoenergetic  and  energy  distributed  sources  with  various 
pitch  angle  dependences. 

Walt  et  alf  were  the  first  to  provide  altitude,  energy,  and  pitch  angle  infor¬ 
mation  on  auroral  electron  fluxes  in  the  KeV  range.  To  do  so  they  obtained  a 
numerical  solution  of  a  Fokkcr-Planck  equation  which  was  originally  used  to  study 

7 

the  properties  of  electrons  trapped  in  the  radiation  belts.  The  method  assumes 

8 

continuous  energy  loss  and  small  angle  scattering.  Banks  et  al,  joined  together 

Walt's  Fokker- Planck  method  and  a  low  energy  approx  ini  ate  two -stream  transport 

9 

method  by  Banks  and  Nagy,  "  previously  applied  to  photoelect ron  tr  ansport.  The 
resulting  equation  was  solved  numerically. 

Beger  et  al,  chose  to  examine  auroral  electron  transport  by  applying  Monte 
Carlo  techniques.  ^  Backscatter  yields,  backscatter  spectra,  and  altitude 
profiles  of  the  energy  deposition  rate  are  among  the  transport  quantities  that  were 
calculated  by  this  method.  In  the  first  of  the  two  papers  noted  above,  information 
was  also  given  on  the  lateral  spreading  of  KeV,  narrow  electron  beams. 

2.  Spencer,  L.  V.  (1955)  Theory  of  electron  penetration,  Phys,  Rev.  98:1597. 

3.  Spencer,  L.  V.  (1959)  Knergy  dissipation  by  fast  electrons.  National  Bureau 

of  Standards  Monograph  1 .  ~ 

4.  Rees,  M.  H.  (1983)  Auroral  ionization  and  excitation  of  incident  energetic 

electrons,  Planet,  Space  Sci.  11:1209. 

5.  Grun,  A.E.  (1957)  Lumineszenz-photometrische  irm^sungen  der  energic- 

absorption  im  st  rahlungsfeld  von  elect  ronquellen,  cmdimensionaler  fall  in 
lu ft,  Z.  Naturforsch. ,  Ser.  A,  12:89. 

8.  Walt,  Al.,  MacDonald,  W.  M. ,  and  Francis,  W.  E.  (1987)  Penetration  of 
auroral  electrons  into  the  atmosphere.  Physics  of  the  Alagnetosphere, 

R.L.  Carovillano,  J.F.  McClay,  and  H.  R.  Radoski,  eds. ,  D.  Reidol, 
Dordrecht,  Netherlands,  534. 

7.  AlacDonald,  W.  Al.  and  Walt,  Al.  (1981)  Distribution  function  of  magnetically 

confined  electrons  in  a  scattering  atmosphere,  Ann,  Phys.  1 5:44. 

8.  Banks,  P.  Al.,  Chappell,  C.R.,  and  Nagy,  A.F.  (1974)  A  new  model  for  the 

interaction  of  auroral  electrons  w’ith  the  atmosphere:  Spectral  degradation 
backscatter,  optical  emission,  and  ionization,  J.  Geophys.  Res.  79: 1459. 

9.  Banks,  P.  Al.  and  Nagy,  A.F.  (1970)  Concerning  the  influence  of  elastic 

scattering  upon  photoelectron  transport  and  escape,  J.  Geophys.  Res. 
75:1902.  ~ 

10.  Berger,  Al.J.,  Seltzer,  S.  Al.  ,  and  Alaeda,  K.  (1970)  Energy  deposition  by 
auroral  electrons  in  the  atmosphere,  J.  Atmos.  Terr.  Phys.  32:1015. 

11.  Berger,  Al.J.,  Seltzer,  S.A1.,  and  Alaeda,  K.  (1974)  Some  new  results  on 
electron  transport  in  the  atmosphere,  J.  Atmos.  Terr.  Phys.  38:591. 


There  are  at  least  throe  transport  models  utilizing  the  linear  transport 

equation  which  have  recently  been  applied  o  auroral  studies.  These  are  the 

l ?  i  *■}  14 

models  of  Strickland  et  al,  “  Mantas  and  Stanmes.  All  three  models  give  a 
detailed  description  of  elastic  scattering  and  alio  a-  for  discrete  energy  loss.  Dif¬ 
ferences  arise  in  the  representation  of  the  flux  within  the  collision  integral  and  in 

the  method  of  integration  over  depth.  Strickland  et  al,  allow  the  flux  to  vary 

12 

quadrat  ically  in  {  n  K  and  .i nearly  in  y  within  any  given  F,  p  cell.  “  The  quadratic 
dependence  was  introduced  because  of  energy  conservation  problems  for  a  linear 
dependence  when  treating  energetic*  fluxes  above  several  KoV.  The  integration 
over  depth  is  carried  ou‘  by  either  a  finite  difference  method  (second  order 
pred 'etor-cor rector)  or  an  eigenvalue  method.  The  latter  approach  was  found  to 
be  mu  'h  faster  and  more  accurate.  Mantas  allows  the  flux  to  vary  linearly  in 
both  K  and  p  within  any  K,  p  cell.  He  also  uses  a  linear  dependence  within  a 
given  z  cell  which  leads  to  the  standard  finite  difference  expression  for  the  first 
order  z  derivatives.  Stamnes  treats  the  p  dependence  of  the  problem  by  the  dis¬ 
crete  ordina*e  method  and  considers  the  flux  to  be  constant  within  a  given  F  cell 

14 

(commonly  called  the  multigroup  approximation).  For  a  good  discuss  ion  <if  the 

lr>  l‘> 

discrete  ordinate  and  multigroup  methods  see  Davison.  '  Like  Strickland  et  al, 

14 

Stamnes  uses  an  eigenvalue  technique  to  carry  out  the  integration  over  depth. 

The  discrete  ordinate  method  for  treating  the  angular  dependence  of  the  electron 
flux  converges  poorly  for  highly  anisotropic  scattering  ker  nels.  Therefore,  the 
angular  dependence  of  the  results  of  Stamnes  at  high  energies  in  the  backward 
direction  are  questionable.  All  of  the  transport  methods  lead  to  a  truncated 
matrix  equation  which  is  solved  numerically. 

1  nc  new  feature  of  the  work  presented  by  us  in  this  paper  is  that  we  give 
approximate  analytic  solutions  for  the  primary  auroral  electron  flux  in  the  down¬ 
ward  hemisphere  and  related  quantities.  All  previous  work  on  this  problem  has 
resulted  in  numerical  solutions.  The  analytic  results,  though  approximate,  are 
useful  in  giving  insight  into  the  physics  of  auroral  electron  precipitation  and  in 
carrying  nut  further  analytical  studies  of  auroral  phenomena.  For  example,  in 
work  currently  in  progress  using  the  multiple  scattering  method,  we  show  that 
elastic  scattering  can  be  included  rigorously  and  that  the  solution  given  in  this 


12.  Strickland,  D.J.,  Book,  D.  L. ,  Coffey,  T.P.,  and  Fodder,  J.A.  (197(>) 

Transport  equation  techniques  for  the  deposition  of  auroral  electrons, 

J.  Goophys.  Hes.  81:2755. 

13.  Mantas,  G.  P.  (1975)  Theory  of  photoelectron  thermalization  and  transport 

in  the  ionosphere,  Planet.  Space  Sci.  23:337. 

14.  Stamnes,  K.  (1978)  A  theoretical  investigation  of  the  interaction  of  auroral 

electrons  with  the  atmosphere,  Ph.  D.  thesis,  lTni versify  of  Colorado, 
Colorado. 

15.  Davison,  B.  (1957)  Neutron  Transport  Theory,  Oxford  Press,  London. 
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paper  is  just  the  leading  term  in  an  iterative  sequence  of  solutions  which  treats 
the  complete  problem  including  secondary  electron  production.  In  work  also  in 
progress  the  method  presented  in  this  paper  is  used  to  solve  the  coupled  proton- 
hydrogen  precipitation  problem  for  the  proton  aurora.  These  Solutions  for  the 
proton  and  hydrogen  flux  and  related  quantities  a  re  similar  to  the  ones  presented 
here  for  the  election  precipitation  problem. 

2.  BASIC  KQl.  ATI  ON 

For  a  one  constituent  atmosphere  the  linear  transport  equation  for  energetic 
electrons  in  plane-parallel  geometry  is 


u  *  n(z)Q  (K)  $(/,  K,  y)  n(/)Q 


t<i:>  J  ciK1  (!«'  nu:v  •  i:«)  $(/, 


where  cf>  is  the  electron  flux,  n  is  the  neutral  density,  and 


K“)s  qTk)  I  0i(KV  •  **•>  ■ 


Here,  /  is  the  altitude,  K  the  particle  energy,  /u  the  cosine  of  the  angle  between 
the  particle  velocity  and  the  positive  /.-axis,  and  cn  the  differential  cross  section 
per  unit  energy  range  (referred  to  hereafter  as  the  differential  cross  section)  for 
electron-neutral  particle  scattering  where  the  neutral  particle  makes  a  transition 
from  the  ground  state  to  the  final  state  i.  The  total  c'ross  section,  Q  (K),  is  re¬ 
lated  to  the  differential  cross  section  bv  the  formula 


)t(K)  =  2a  j  ciK'  (V  X  a.(Ka  -  KV)  . 


:t.  K(.M  \TION  KOR  TIIK  PRIM AIO  M.KC.TRON  H.l  \ 


The  differential  cross  sections  describe  elastic,  excitation,  and  inm/ation- 
type  collisions  between  the  electrons  and  the  neutral  particles.  In  ioni/at  ion-type 
collisions  additional  (secondary)  electrons  are  produced.  The  electron  flux  in 
Kq.  (1)  includes  both  primary  and  secondary  electrons.  The  secondary  electron 


spectrum  is  highly  peaked  at  low  energies  (<K  -n)  with  very  few  secondary  elec¬ 
trons  at  higher  energies  *  In  the  auroral  ionosphere  Kmm  somewhat 

arbitrary  but  is  on  the  order  of  1  2  KeV  for  primary  elet  Irons  ranging  from 
1  to  hO  KeV.  If  we  restrict  the  energy  range  of  Kq.  (1)  to  high  energies,  we  may 
approximate  the  electron  flux  by  neglecting  the  secondaries  compart'd  to  the 
primaries;  that  is, 


primaries 


k  ^  i: 


(4) 


and  use  a  differential  cross  section  for  ioni/at  ion-type  collisions  that  deset  ibc-. 
only  the  scattering  of  the  primary  electrons.  Viewing  the  ioni/ation-tvpe  colli¬ 
sions  in  this  wav  we  may  approximate  the  differential  cross  sections  for  KeV 
auroral  electrons  by 


Ell)  C2*)~LQc(E'n(E'  -  E) 


(a) 


i:tl)  '  <2a-»-1  cioxk<i-:*) c  -  (i:  -  w,.)l  p,>xk<u\u' 


(Hi 


'”<i:v  -  Ki,)  '  (2»C1a|P>  O  ’.  O  pjfhw'.u)  . 


i  ( 


(7) 


\\  , 


Merc,  the  subscripts  c,  c\,  and  i  refer  to  elastic,  excitation,  and  ioni /at  ion -t  ypi 
collisions,  respecM  ively,  6  is  the  Dirac  delta  function,  the  p  functions  are  sym¬ 
metric  m  u 1  and  u  and  are  normalized  to  one,  t lit'  total  cross  section  for  each 
process  is  denoted  by  Q,  and  tilt*  superscript  p  refers  to  primary  elect  rolls 
is  the  excitation  energy,  Kor  a  discussion  f  '!,*  'dnstio,  excitation,  and  rr'i. 
t  ion  type  cross  sections  sec  Strickland  et 

The  equation  for  the  primary  electron  flux  is  Kq.  (1)  with  the  d  i  ffe  rent  i  a  1 

f71  and  is  valid  for  1!  >  K 

m  m 


sections  given  by  Kqs.  (a),  (h),  and  IT),  and  is  valid  f » w  1.  >  K 


\  KOKW  MU)  SC  \TTKK1M.  \M>  W  KR  \CI  DISCRKTK  F.M  H(i) 

I.OSS  \PPRO\IM  VTIONS 

In  the  auroral  region  two  facts  about  the  differential  cross  sect  ions  at  high 
energies  are  apparent.  Thev  are  that  1)  the  differential  itoss  sect  ions  ;nv  highl\ 
peaked  in  the  forward  direction,  the  excitation,  and  iom/ation  tvpe  collisions  ben 


DC  (ipal,  (  .IK.  I’eterson.  W.  K.  ,  ami  Heaty.  K.C  .  (lbTU  Measu  for  lent  s  of 
sec<  inda  rv -elect  ron  spectra  produced  b\  elect  ron  impact  tom  rat  :•  »r,  . 

number  of  simple  gases,  J.  Cliem.  1’hys.  c>:4100. 


even  more  highly  peaked  than  the  elastic  collisions,  and  2)  the  average  energy 
loss  per  inelastic  collision  (excitation  and  ionization)  is  a  weak  function  of  the  in 
cident  electron  energy.  This  suggests  that  for  inelastic  collisions  we  make  the 
forward-scattering  and  average  discrete  energy-loss  approximations: 


k  \k 


exkv 


(E«>\  <5lE»  -  (E  +  W)1  60i‘  - 


(8) 


£  -  En>  -  (2*)'1  (Y  (E')\  6[E'  -  (E  +  W)J6(m'  -»)  ,  (9) 


t 


where  W  is  the  average  primary  electron  energy  lost  per  inelastic  collision 
defined  by 


W  h 


r  i  wk<wE> +  z  /dE,(K  - 


-l 


‘I  Qexk<E)+  Z/dE’a!/>(E*E,)‘ 

-  k  f 


(10) 


W  as  a  function  of  E  for  the  auroral  cross  sections  used  by  Strickland  et  al,  is 


32,  36,  42  and  48  in  eV  for  1,  3,  10  and  30  KeV  primary  electron  energies, 

12 

respect ivelv.  In  the  average  discrete  energv-loss  approximation  \Y  is  assumed 


constant.  For  elastic  cross  sections  a  screened  Rutherford  cross  section  is 
often  used: 


oe(  EV  -  Eu)  ^  (2^)*’iQe(E')6(E«  -  E>  pfa,  *i\ *i>  , 


(11) 


pin,  w',p) 


2 p( l  +_r) Hi _+  2 n  -  /j y ) 


Kl  +  2n  -  -  (l  -  »'2)(\  -  u 2)13/2 


(12) 


where  r)  is  a  function  of  energy.  For  graphs  of  p  as  a  funtion  of  q,  (jt  and  ^ ,  and 

1 2 

n  as  a  function  of  E  see  Strickland  et  al. 


Inserting  Kq.  (8)  through  Kq.  (12)  into  Kq.  (1)  and  t  ransf<  mining  Kq.  (1)  t< 
an  equation  m  terms  of  the  optical  depth,  t,  where 

dr  ■*  -n(/.)Q(K)  dz  ( 

we  obtain 

(-U  »  lj  $ ><T.  k.m>  ^-Q{-yy--  i:  *•  W),  i:  ••  u.  M) 


Q  ,<K) 


-  4>(T,  -  J  dw'  <1><t,  K,mM 


where  Q  and  denote  the  total  inelastic-  and  total  elastic  cross  sections,  respee 
tively,  and  t  denotes  t (z,  K).  The  ratio  Q(K‘  *  W)  ’Q(F)  is  near  1  and  the  ratio 

Q  (K)  Q(K),  as  given  by  Strickland  et  ai,  is  0.84,  0.  82,  0.45  and  0.27  for  1,  2, 

°  1  '> 

10  and  20  KeV  primary  electron  energies,  respectively.  “  On  the  right-hand  sid 

of  Kq.  (14)  we  may  view  the  second  and  third  terms  as  producing  a  correction  to 
tlie  first  term.  This  cor  rection  is  not  necessarily  small.  Fur  example,  we  know 
that  in  the  backward  hemisphere  elastic  sc  attering  is  the  dominate  process  and 
contributes  the  most  to  the  backscattered  flux.  The  solution  here  is  very  sensi¬ 
tive'  to  the  ratio  Q  'Q  and  to  the  precise  shape  of  p(n, iu m).  However,  in  the' 
downward  hemisphere  we  expect  the  first  term  on  the  right-hand  side  of  Kq .  (14) 
to  product'  a  contribution  to  4>  which  is  comparable  to  or  larger  than  the  contribu¬ 
tion  made  by  the  second  and  third  terms.  In  this  paper  we  seek  an  approximation 
to  the  flux  in  the'  downward  hemisphere.  An  equation  for  such  an  approximate 
solution  is  obtained  from  Kq.  (14)  by  taking  the  limit  of  the  function  p  as  K  -*  oc 
(q  -  0).  It  can  be  shown  that 


pin,  v  \  ijl)  -+  6((jl  -  (j) 

n  +  0 


In  the  high-energy  approximation  elastic  scattering  becomes  highly  peaked  in 
the  forward  direction  and  drops  out  of  Kq.  (14).  We  obtain 


;  i,  •  .i 


*(T<z,  E),  K,m)  <f(r(z,  K  <  VV),  K  +  W ,  u)  .  (lii) 


Equation  (16)  is  a  partial  differential-difference  equation  for  the  primary  electron 
flux  now  denoted  by  <t>,  as  a  function  of  t,  K  and/j,  approximately  valid  for 
0  <  T  <  ge,  K  <  K  <  or,  and  -1  <  /j  <  0. 


5.  MILTIP1.K  SCATTERING  MKT  HOI) 

The  multiple  scattering  method  may  be  applied  in  principle  to  any  linear 
transport  equation  satisfying  certain  general  requirements.  For  a  discussion  of 

the  multiple  scattering  method  applied  to  elastic  scattering  see  Goudsmit  and 

17  IB  19 

Saunderson  and  Wang  and  Guth,  to  inelastic  scattering  see  Kano,  and  to 

20 

both  types  of  scattering  see  Case  and  Zweifel.  “  In  our  case  we  wish  to  apply  it 
to  Kq.  (16)  which  contains  only  inelastic  scattering.  To  do  so  we  rewrite' 

Kq.  (16)  as 


-m  §7  *  l)  * 


& 


(17) 


where  d[  is  the  linear  operator  that  shifts  the  energy  variable  of  <!>  by  W  and 
multiplies  bv  Q(K  K  W)/Q(K).  The  multiple  scattering  method  consists  of  writing 


<t(T,  ft,  n )  -  Y,  4>n<T.  I'.m) 

n-0 


(18) 


where  the  satisfy  the  infinite  system  of  equations 


*  W  +  1  *o  °  ' 


(19) 


3T 


+  1  j  <t>  l  £<b0  , 


(20) 


17.  Goudsmit,  S.  and  Saunderson,  J.  L.  ( 1940)  Multiple  scatter ing  of  elect rons, 

Phys.  Rev.  57:24. 

18.  Wang,  M.C.  and  Guth,  E.  (1951)  On  the  theory  of  multiple  scattering,  par¬ 

ticularly  of  changed  particles,  Phys.  Key.  84:1092. 

19.  Fano,  U.  (1953)  Degradation  and  range  straggling  of  high-energy  radiations, 

Phys.  Rev.  92:330. 

20.  Case,  K.  M.  and  Zweifel,  P.  L.  (1967)  Linear  Transport  Theory,  Addison- 

Wesley,  Heading,  48.  ““  ~~ 
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{21) 


pn- 1 


and  then  solving  for  the  4^  subject  to  appropriate  boundary  conditions.  If  $  exists 
it  is  referred  to  as  the  multiple  scattering  or  orders  of  scattering  solution.  The 
fact  that  Eq.  (18)  is  a  solution  to  Eq.  (17)  is  seen  by  adding  the  above  infinite 
system  of  equations  for  the 

The  function  4^  which  we  may  call  the  nth  order  partial  flux  has  a  simple 
physical  interpretation.  is  the  flux  of  particles  per  unit  energy  per  unit  solid 
angle  at  E,  and  ^  which  have  experienced  no  collisions,  ^  is  the  flux  of 
particles  per  unit  energy  per  unit  solid  angle  at  z,  E,  and  ju  which  have  scattered 
once,  and  so  on  for  each  order  of  scattering  denoted  by  n. 


6.  MULTIPLE  SCATTERING  SOLUTION 


The  boundary  condition  for  electrons  incident  at  the  top  of  the  ionosphere  is 
$  (0,  E,  =  $  (E,  p)  ,  for  - 1  <  ju  <  0  (22) 

This  implies  that  for  n  =  0, 

4  (0,  E,»<)  »  «  (E,p)  ,  for  "lSfi  <  0  ,  (23) 

and  for  n  >  1, 


</>n(0,  E,^)  =  0  ,  for  -1  <  m  <  0 


(24) 


No  particles  are  incident  from  below  the  ionosphere. 

The  solution  to  Eq.  (16)  may  be  found  as  follows.  The  solution  for  <j>  subject 
to  the  above  boundary  condition  is 


$(E,/u)  exp  (t/m)  ,  -1<  fjL  <  0 


4q(t,  E,u)  = 


,  0  ^  1 


(25) 
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and  the  solution  <£>n  in  terms  of  for  n  1  is 

T 

<^n(T,  E#Ju)  =  -  exp  (t/^)  J  6 1  /lx  1  bi  i  K  +  exP  <  “t  /m)  ,  (2*»> 

o 

for  - 1  <  jli  <0.  For  0  <:  fu  <  1  the  art*  zero.  The  complete  solution  is 

00 

$(t,  !•:,  M)  =  <J>(E,  a)  exp  (t/m)  +  £  $>(K  +  nW,  /j)  b^E) 

n- 1 

x  [nn(T/n,  i:)  -  nn(o,  »•:>  exp  f  (27) 

for  -1  <  ju  0,  and  $  =  0  for  0  ^  /a  <  1.  The  recursion  relation  for  the  H  func¬ 
tions  is 

exp  <-x)  Hn+1(x,  K)  -  Hn+1(0,  K)  = 

X 

-  f  dx'  exp  (-x')[Hn(b11x>,  E  +  W)  -  IIn<0,  K  +  W)  exp  (bnx')l  .  (28) 

o 

An  explicit  formula  for  t lie*  general  H  function  ha s  been  worked  out  bv  mathemat¬ 
ical  induction  but,  as  it  is  quite  lengthy,  only  the  first  three  functions  will  be 
given  here: 


Hjfx,  1C) 

iijj  exp  (bnx)  . 

(29) 

11., (x,  1C) 

a22a21  oxp  (b>lx)  *  a22:,ll  oxp  <bUx)  • 

C10> 

H.<(x.  1C) 

5  a33a33a3I  °XP  (b:Ux)  ~  a33a32al  1  oxp  <bllx> 

*  a3.ja22  '21  °xp  (b21x)  +  a33a22a  1 1  exp  (bllx)  ■ 

cm 

the  above 

forimilao  we  have  introduced  the  following  definitions: 

1)  h 

mu 

nnl(K)  -  Q(IC  +  n\V  )  'Q(  1C  +  (ni  -  1>W>  , 

C12) 

Hi 


Mill 


Equation  (27)  may  be  rewritten  in  an  alternative  form  as 


oo  n 

$(t,  E,u)  =  $>(E,*j)  exp  (t/jj)  +  $(E  +  nW,^)  £  A 

n=  1  1  =  1 


X  [exp  (Bj  t/u)  -  exp  (t/^)]  , 

(34) 

where 

B^E)  =  Bi  s  b^E)  , 

(3  5) 

and  the  first  few  =  A^(E)  are 

An  =  +  bnan  - 

(3f>) 

A22  =  +  b2 la22a2 1  ■ 

(37) 

A21  =  "  b2 la22al 1  * 

(38) 

A33  =  +  b31a33a32a31  ’ 

(39) 

A32  =  "  b3 la33a22a2 1  ’ 

(40) 

A3  1  =  '  b31  a33a32a  1 1  "  a33a22all^ 

(41) 

Several  comments  about  this  solution  are  in  order.  The  primary  electron 
flux,  is  an  exact  solution  to  Eq.  (16)  which  is  in  turn  an  approximation  to 
Eq.  (1).  The  fact  that  it  solves  Eq.  (1(5)  exactly  can  be  verified  bv  direct  sub¬ 
stitution.  The  value  of  $  depends  linearly  on  its  boundary  value;  doubling  the 
incident  flux  doubles  the  flux  at  all  altitudes.  The  solution  yields  no  backscattered 
primary  electron  flux.  The  reason  for  this  is  that  in  the  forward  scattering  ap¬ 
proximation  elastic  scattering  drops  out  and  no  particles  are  scattered  into  the 
backward  direction. 


r 


7.  HEMISPHERIC  ALLY  AVERAGED  FLUX,  ENERGY  DEPOSITION 
RATE  AND  IONIZATION  RATE 

The  hemispherical  electron  flux,  is  defined  as  the  average  value  of  the 

electron  flux  over  the  downward  hemisphere  and  is 

o 

$h<t,  E)  =  2?r  J  d*j  *(t,  E,ui)/2n  .»  (42) 

-1 

The  energy  deposition  rate,  is  defined  as  the  energy  deposited  per  unit  vol¬ 

ume  per  unit  time  by  the  precipitating  electrons,  and  is  given  in  plane-parallel 
geometry  by 


oo  -+T 

nE(z)  s  2-n  J  dE  J  d„MKn(z)Q(E)  4>(t,  E,*i>  ■  (43) 

o  - 1 

The  ionization  rate,  n^,  is  defined  as  the  number  of  electron-ion  pairs  produced 

per  unit  volume  per  unit  time.  This  can  be  found  to  a  good  approximation  from 

2 1 

bv  using  Bethe's  formula 

h, 

Ojb.)  nK(*)/Ks  ,  (44) 

where  E  is  a  constant,  A  good  value  for  E  for  auroral  electrons  has  been 
s  s 

found  to  be  34  eV.  In  this  way  we  find  an  approximate  formula  for  rj.  without 
solving  explicitly  for  the  secondary  electron  flux. 


8.  SOLI  TIONS  FOR  SPECIFIC  INCIDENT  ELECTRON  FLUXES 
8.1  Preliminary  Remark* 

In  this  section  we  give  analytic  expressions  for  the  quantities  defined  in 

Section  7  for  several  specific  forms  of  the  incident  electron  flux*  Before  we  do 

so,  we  remind  the  reader  that  our  analytic  solution  for  the  electron  flux  is  only 

valid  for  E  >  E  .  *  However,  in  the  energy  integral  from  0  to  oc  which  defines 
m  m  * 

r}^,  we  note  that  the  contribution  to  the  integral  below  E  only  amounts  to  a  few 
'E'  ®  mm 

percent  for  auroral  electrons. 


21.  Bethe,  H.A.  (1933)  Iianbuch  dor  Physik,  Verlag  Julius  Springer,  Berlin 
24:491. 


8.2  U nidi reetional-Moiioener^e tic  Incident  Flux 

For  a  unidirectional -monoonorgotio  electron  flux  incident  at  the  top  of  the 
ionosphere  the  boundary  condition  is 


$  (E,  p ) 


6<E  -  Eq)6(m  +po) 


(4f>) 


where  Qg  is  the  total  energy  flux  in  the  downward  direction,  E^  is  the  electron 
energy,  and  is  the  cosine  of  the  angle  of  incidence,  x,  where  05  x  <  ff/ 2. 
The  electron  flux  is 


$(t,  E,m) 


Q,. 


2  ttE  u. 

oMo 


E^)  exp  (t/jh) 


+  y  6(E  +  nW  -  E  )b  .(E)  [H  (r/M,  E)  -  H  (0,  E)  exp  (t  /^)]  >  . 

o  n  i  n  n  j 

n=l 


The  hemispherical  flux  is 


(48) 


Q 


*h(t’k)=  [2^;Ji6(E  -  v exp  ('T/wo) 


+  y  6(E  -  (Eq  -  nW))bnl  fiy-T/^,  E)  -  Hn(0,  E)  exp  <-T/no)l) 


n=l 


The  energy  deposition  rate  is 


nE<^=  Q<Vn(z> 


jexp  (-t(z,  E 0>/u0>  +  Y,  ^  '  nW/^o>  b^j(Eo  -  n\V) 
^  n  =  l 


X  [I  (-r(7.  E  -  nW )/u  ,  E  -  nW)  -  H  (0,  E  -  nW) 

n  o  ~  ~ 


X  exp  ( “T(z,  -  nW)/*i^)]j> 


(48) 
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where  the  I  functions  are 


In<*-E>*  SxHn<*'E>  • 


It  should  be  pointed  out  here  that  the  above  sums  do  not  extend  to  oo  but  to  a  max¬ 
imum  value,  N,  given  by  NW  =  E  -  E'  ;  . 

'  &  J  o  min 


8.3  Isotropic- Maxwellian  Incident  Flux 

For  an  isotropic -Maxwellian  electron  flux  incident  at  the  top  of  the  ionosphere 
the  boundary  condition  is 


$(E,p)  =  - E  exp  (-E/E  )  ,  for  -1  <  /u  <  0 


where  Q  is  the  total  energy  flux  in  the  downward  direction  and  E  is  the  charar- 
s  o 

teristic  energy.  The  resulting  electron  flux  is  given  by 


$(t,  E,/u)  =  Ve  exp  (-E/Eq)  exp  (t/^)  +  £  <E  +  nW) 


exp  (-  (E  +  nW)/Eo)  bnl  [H n<r/M,  E)  -  Hn<0,  E)  exp  {t/u)\ 


with  $  =0  for  0  <  ^  "S  1,  The  hemispherical  flux  is 


$ „(z,  E)  =  I  ^ 
H  \2?r  o 


f  00 

<E  exp  (-E/Eo)  E2(t)  +  £  (E  +  nW) 


exp  (-(E  +  nW)/E  )  b  ,  [K  (t,  E)  -  H  (0,  E)  E„(t)1  \  , 

'on  in  n  2  J 


where  E^Cx)  is  the  second  exponential  integral  function  defined  by 


E9(x)  J  dt  t  c  exp  (-xt) 
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The  first  three  K- functions  are 


K^x,  E)  a11E2(bnx)  ,  (54) 

K2(x,  E)  =  a92a21i:2(b21x)  -  a99a1  jE2^bi  lx^  *  (55) 

K3(x,  E)  “  a33a32a31L2^b31X^  "  a33a32al  1  ^2^1 1X* 

’  a33a22a21^2^b21x^  +  a33a22a 1 1 1  J2(b 1 1X*  *  ^5(>* 

The  general  rule  for  the  K  functions  is  the  same  as  that  for  the  H  functions  with 
the  exponential  function  replaced  by  the  second  exponential  integral  function. 

The  energy  deposition  rate  is 


”e(z) 


oo 

n(z)  f  dEEQ(E) 
E 

min 


E  exp  (-E/E  )  E2(t) 


oo 

+  £  (E  +  nW)  exp  (-(E  +  nW)/E0)bnllLn(T,  E)  -  Hn(0.  E)E  (T)} 


The  rule  for  the  L  functions  is  the  same  as  that  for  the  I  functions  except  that  the 

exponential  function  is  replaced  by  the  second  exponential  integral  function.  In 

performing  the  integral  over  the  total  inelastic  cross  section,  Q,  the  integral 

mav  be  truncated  at  some  E  .  When  this  is  done  the  infinite  sum  terminates 

max 

at  a  maximum  value,  N,  given  by  NW  =  E  -  E  .  . 

&  J  max  min 

8.4  l«otropic-Monoenergeti«  Incident  Flux 

For  an  incident  isotropic -monoenergetic  flux  the  boundary  condition  is 


$(E, u) 


-1  <  m  <  0 


(58) 


where  Q  is  the  total  energy  flux  in  the  downward  direction  and  E  is  the  electron 
®  o 

energy.  The  expression  for  is 
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$(t,  E,/u)  = 


“'I 


E  )  exp  (t/ih)  + 


I 

n-1 


6(E  +  nW  -  K  )b  , 
o  nl 


X  [Hn(T//,,  E)  -  I!n(0,  E)  exp  (t/m)]  j  ,  (59) 

with  $  0  for  0<  M  1. 

9.  RELATIONSHIP  RETAKEN  THE  ENERGY  DEPOSITION  RATE 

AND  CHAPMAN  AND  J  El  NOTIONS 

9.1  Prcliiuiiiar)  Remark# 

Chapman  functions  are  functions  of  altitude  which  arise  in  the  Chapman  theory 
of  the  daytime  ionosphere.  The  functions  gives  the  ionization  rate,  which  is  re¬ 
lated  to  the  energy  deposition  rate,  for  a  single  Chapman  layer.  The  ionization 
rate  profile  is  produced  by  (he  photoionization  of  the  neutral  particle  gas  by 
ionizing  elect  romagnet  ic  radiat  ion  from  the  sun.  In  the  auroral  zone  at  night, 
the  ionized  layer  is  produced  by  a  different  process;  the  precipitation  of  energetic 
electrons  incident  at  the  top  of  the  atmosphere.  It  has  been  known  for  some  time 
that  the  auroral  ionization  rate  profile  (or  energy  deposition  rate  profile)  is  much 
more  highly  peaked  as  a  function  of  altitude  than  is  a  single  Chapman  function. 

In  this  section  we  establish  the  connection  between  the  Chapman  functions  and 
the  auroral  ionization  rate,  specifically  we  show  that  the  solution  for  the  energy 
deposition  rate  in  the  auroral  ionosphere  is  a  superposition  of  Chapman  functions 
for  unidirectional  incident  electron  fluxes,  and  a  superposition  of  what  we  call 
J  functions  for’  isotropic  incident  electron  fluxes.  Tn  order  to  do  this  wo  use  the 
alternative  expression  for  rj ^ .  we  obtain  when  Eq.  (34)  is  substituted  into  Eq.  (43). 

00  o 

nF(z)  2*  n(z)  J  dE  J  d*i  EQ(K)  j$(E,/u)  exp  (r(z,  1-:)  u) 

E  .  - 1  ^ 

m  in 

oo  n 

+  y  $o:  +  n\v,n)  y 

n  1  l  - 1 


V  [Bf(K)  exp  (Bf(E)T(z,  K)/ju) 


exp  (t(z,  E) 


<<>0) 


V 


9.2  Generalized  Chapman  and  J  Functions 


The  Chapman  functions  as  introduced  by  Chapman  is 
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Ch(z)  =  n(z)  exp 


(61) 


where  f  is  a  constant  and  ^  is  the  cosine  of  the  angle  of  incidence,  x»  where 
0  s  x  <  */2.  Let  us  now  define  a  generalized  Chapman  function  of  three  argu¬ 
ments 


n(z)  exp 


HE) 

JU 


00 

J  dy  n(y) 

z 


(62) 


where  0  <  z  <  oo,  -1  S  ju  <  0,  and  f  is  a  positive  function  of  E,  where  0  <  E  <  oo. 
By  analogy  with  the  Chapman  function  let  us  introduce  a  J  function 


o 

00 

f  dfu  n(z)  exp 

-  f  dy  n(v) 

1 

z 

and  a  generalized  J  function 


(63) 


o 

oo 

J(z,  f(E))  s  J  d/u  n(z)  exp 

f  dy  n(y) 
z 

-1 

(64) 


where  Os  z  <  oc,  and  f  is  a  positive  function  of  E,  where  Os  E  <  oo.  Transform 
ing  Eq.  (64)  to  an  integral  on  t,  where  t  -  -/u  \  we  see  that 


J(z,  f(E))  -  n(z)  E2 


oo 

f(E)  J  dy  n(v) 

7 


(65) 


where  E9  is  defined  by  Eq.  (53). 


22.  Chapman,  S.  (1931)  The  absorption  and  dissociative  or  ionizing  effect  of 
monochromatic  radiation  in  an  atmosphere  on  a  rotating  Earth,  Pt  oc. 
Phys.  Soc.  43:26. 


9.3  Energy  Deposition  Hate  for  l' nidi  rational  Incident  Fluxes 

For  a  unidirectional  incident  flux  cos  0  <  x  <  ff/2)  the  energy  depo¬ 

sition  rate  given  by  Hq.  0>0)  is 

K 

max 

nF(z)  =  2  t  J  dEKQ(K) 

K 

min 

N  n 

*  X  X  *(E  +  n\V)An<(E) 

n  - 1  l  l 


AK)  Ch  (z,  -  —  ,  B#(K)Q(E) 
*  \  ^  i 


(Hn) 

Thus,  the  exact  solution  of  Kq.  (1(0  yields  the  energy  deposition  rate  as  an  energy 
integral  over  u  finite  superposition  of  generalized  Chapman  functions. 

The  explicit  solution  for  a  unidi rect ional-monoene  rgotic  incident  flux  ;> 


N  n 

+  I  I  <1  -  nW/E0)Bf(Ew  -  n\V)Ani(Ko  -  nW) 
n  - 1  f  - 1 


v 


BA  K  -  n\V)Ch 
f  o 


-  — ,  B,(E  -  n\V)Q(K  -  n\Y)] 
^  i  o  o  / 

°  / 


(Ii7) 


This  explicit  solution  reduces  to  a  finite  superposition  of  generalized  C  hapman 
functions  with  shifted  energy  arguments. 


9.1  KntT£\  Deposition  Kale  tor  Isotropic  Ineiilenl  Maxes 

For  ;m  isotropic'  incident  flux  t lit*  energy  deposition  rut*-  us 
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max 
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Hero,  the  exact  solution  of  Hq.  (lfi)  yields  an  energy  deposition  rate  as  an  energy 
integral  over  a  finite  sum  of  generalized  J  functions. 

For  an  isot ropic-IMonoenerget ic  incident  flux  the'  enerv  integral  can  be  doin' 
explicitly  and  for  an  isot  ropie-Maxwellian  incident  flux  the  energy  integral  can  be 
done  by  quadr  ature. 


10.  I  SF  OF  PSKl  DOF  YKTICFFS  TO  \PPR0\1MATK  TDK  SI  MS 

In  the  formulae  of  Sections  (>  through  9.  W  is  a  number’  on  the  order  of  40  eV 
and  the  primary  auroral  electron  energies  range  from  about  1  to  (iO  KeY.  As  a 
result,  the  sums  in  the  formulas  extend  to  as  many  as  lr>00  terms  representing  as 
many  as  1500  scatterings  before  an  energetic  electron  loses  all  its  energy.  How¬ 
ever,  it  turns  out  that  a  good  approximation  to  these  largo  sums  may  be  obtained 
by  introducing  the  notion  of  pseudoparticles  which,  in  turn,  allows  us  to  take 
many  fewer  terms  in  each  sum.  A  pseudoparticle  is  a  particle  which  as  a  cross 
section  W  'YY*  times  smaller  than  the  real  particle  but  has  an  average  energy  loss 
per*  inelastic  collision  W  /YV  times  greater,  such  that  the  product  of  the  two  re¬ 
nin  ins  the  same.  With  ;his  approximation  a  good  answer  is  obtained  for  the  sums 
given  in  Sections  (>  through  9  with  many  fewer  terms  in  each  sum. 

In  order  to  introduce  pseudoparticles  into  the  equations  of  Sections  t>  through 
9  we  simply  replace  Q  by  and  W  by  W^,  that  is, 

Q(K)  -  Q  (H)  s  <\V  W  )  QUO  ,  019) 

P  P 

\V  -  YV  =  (YV  /YV)  YV  .  (TO) 

P  P 


We  have  examined  the  convergence  of  the  pseudoparticle  method  as  we 
approximate  the  sums  bv  a  smaller  and  smaller  number,  N,  of  pseudoscatterings 


spanning  the  energy  range*  from  E  to  E  .  We*  have  found  that  a  good  answer 
r  f.  f,.  f-.  max  mm  ** 

van  be  obtained  with  a  surpr  isingly  small  number  of  pseudosi  attermgs.  For 
example,  for  an  tsot  ropie -Maxwellian  incident  flux  with  K  ranging  between  1  and 
10  KeV  only  10  to  15  pseudoseat  ter  ings  spanning  the  energy  range  arc  needed  to 
obtain  answers  to  w  ithin  a  few  percent  of  the  full  solution  at  most  altitudes  of 
interest. 

In  Figure  1  we  show  the  hemispherieally  averaged  partial  fluxes,  Ojjn,  for  a 

sequence  of  n  values  for  an  isotropic* -Maxwellian  incident  flux  with  K  2  KeV 

2  ° 
and  =  1  erg;cm“s  for  the  case  where  \  12.  That  is  to  say,  12  pseudoscat¬ 
terings  spanning  the  ener  gy  range'  of  interest  from  E  11  K*'Y  to 

K  .  1  KeY.  In  Figure  2  we*  illustrate*  the*  r  ate  of  convergence  of  $Tr  bv  show - 

ing  the  sequence  of  partial  sums  of  ^|jn  foe  the  same  case.  In  Figure  .‘1  we  show 

the*  partial  ionization  rates  (n.  -  nT,  /34)  for  the'  same  ease.  Note  that  for  n. 

1  in  En  ho 

we  obtain  a  J  function  profile*  and  for  17. ^(n  >  1)  we  obtain  functions  with  a  single 

node,  each  of  which  subtracts  from  r\.  at  high  altitudes  and  adds  to  n.  at  low 
'  u>  H  lo 

altitudes.  In  Figure  4  we  show  the  rate  of  convergence  of  r)-  by  showing  the* 
sequence'  of  partial  sums  of  17.  ^  for  the  above  case. 


4>Hn  <Z,E  =  2  5KeV)  (ei/cm2  s  eV  $fl 


Figure  1.  Hemispherieally  Averaged  Partial  Fluxes  for 
an  Isotropic  Maxwellian  with  E  =2  KeV,  0=1  erg / 
rm  s  and  N  =  12  °  s 


Figure  4.  Sequence  of  Partial  Sums  of  Partial  Ioniza¬ 
tion  Hates  for  an  Isotropic -Maxwellian  as  in  Figure  3 


As  a  general  rule  we  find  that  the  lower  the  altitude  the  more  pseudoscatter¬ 
ings  (the  larger  the  N)  covering  the  range  from  F  to  F  .  are  needed  to  nb- 

max  mm 

tain  a  good  answer.  For  example,  we  find  that  in  the  above  cas^*  the  different  *' 
in  the  energy  deposition  rate  for  N  -  12  and  NT  =  20  is  -0.25  percent  at  180  km. 

+0.  8  percent  at  110  km  and  +1.8  percent  at  104  km. 

11.  COMPARISONS 
11.1  Preliminary  Remarks 

In  this  section,  we  compare  our  analytical  results  to  other  transport  calcula¬ 
tions.  Quantities  of  interest  are  the  electron  flux,  its  homisphorically  averaged 
value,  and  the  energy  deposition  rate.  Selected  comparisons  will  be  made  for 
monoenergetic  and  Maxwellian  sources. 

The  particular  models  we  use  for  comparisons  are  the  Hoes,  Range,  and 
Strickland  models.  A  Range  model  was  specifically  developed  for  this  work  and 
to  follow  will  be  a  brief  description  of  its  contents.  We  start  with  the  loss 
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function,  L(E),  which  one  can  obtain  from  the  Bethe  formula  for  energies  above 
several  hundred  eV  or  from  cross  sections  using 


UE>  =  X 


WkQexk(K>  + 


Y  /dE’(K  -  E'JojyNfc:,  K') 


(<  1) 


L(E)  has  been  calculated  from  the  cross  sections  in  Strickland  el  al,  for  \,)t  and 
from  current  tabulations  for  0^  and  0.  The  formula  needed  to  specify  the 
energy  loss  is 


7 

o 

:  Ko-  I  /  ,  (72) 

j  7 

where  is  the  cosine  of  incident  angle,  z  is  aaitude,  E  is  the  starting  ener  gy 
at  z  ,  and  n<(z)  is  th€»  partic  le  density  of  the  jth  neutral  species.  The  energy 
deposition  rate  is 

p.,(z)  =  J  d  E  J  dn  7  n.(z)$(E  ,/u  )  L.(E(z,  E  ,  fj  ))  .  (73) 

h  •  o  •  o  ^  j  o*  o  i  *  o  Mo 

j 

For  all  four  models  compared  in  this  section  the  Jnrehia  (1000°)w^  model 
atmosphere  was  used.  For  the  Strickland  and  Range  models  the  same  individual 
inelastic  cross  sections  and  constituent  neutral  densities  were  used  and  for  the 
analytic  model  the  same  total  inelastic  cross  section  was  used.  The  weighted 
total  inelastic  cross  section,  Q(K),  was  2.83,  1.31,  0.  0.221,  0.0893,  and 

0.0297  times  10  ^  in  cm"  for  0.3,  l,  3,  10,  30,  and  100  KeV,  respectively. 

The  value  used  for  Ey  was  34  eV  and  the  values  used  for-  \V  are  given  in  Section  4. 

11.2  C^oniparisona  for  Isotropic -Monocncrgotic  Incident  Fluxes 

The  models  to  be  considered  are  the  analytic,  Hees,  and  Range  models. 
Incident  energies  to  be  considered  are  2,  5,  and  10  KeV.  Figure  5  shows  vs  / 
for  incident  power  densities  of  1  erg/cni“s.  The  Rees  results  were  generated  by 
us  using  his  energy  dissipation  function.  The  analytic  and  Range  results  are  sim¬ 
ilar  as  expected  since  the  cor  responding  models  both  contain  the  forward 


23.  Jacchia,  L.G.  (1977)  Thermospheric  temperature,  density,  and  compos. i  ion: 
new  models.  Rep.  375,  Smithson.  Astrophvs.  Observ. ,  C  ambridge,  Mass. 
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10  Kev 


Figure  5.  Energy  Deposition  Hates  from  the 
Analytic,  Range,  and  Rees  Models  for  10,  f>, 
and  2  KeV  Isotropic -Mono^nergetic  Sources 
Each  Containing  1  erg/cm^s 
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scattering  approximation.  Departures  result  from  the  differences  in  the  descrip 
tion  of  energy  loss.  The  analytic  model  allows  f<>  discrete  energy  loss  and  in 
turn  straggling.  Greater  differences  exist  between  the  analytic  and  Rees  results 
the  latter  showing  a  broaden*  deposition  prof''  with  less  total  energy  deposition. 
We  expect  these  latter  results  to  be  the  more  accurate*  since*  they  include*  the 
effect  of  pitch  angle  scattering. 

Il.il  (ioniffctrisotis  tor  Iso  tropic- -Maxwellian  Incident  Kluxes 

The  models  to  he  considered  are  the  analytic.  Range,  and  Strickland  models 
Characteristic  energies  to  be  considered  are  1,  2,  and  5  KeV.  Figure  f>  shows 

the  ine'ide^nt  isotropic-Aluxwcllian  distribut  ions  for  these*  energies.  Again,  the* 

2 

incident  power  densities  are  1  erg/cnTs,  Figure  7  shows  vs  */  for  the  three* 
models  and  three  incident  fluxes.  These  results  exhibit  basically  the  same 
behavior  as  those  appearing  in  Figure*  5.  The  Strickland  results  give*  a  broader 
distribution  with  a  smaller  total  energy  deposiion  rate  due  to  pitch  angle  scatter 
ing,  contained  in  his  model.  For  all  the  analytical  results  of  this  section,  12 
pseudoscatterings  were  used  spanning  the  range  from  1  to  13  KeV.  Based  upon 
our  studies  12  pseudoscatterings  give  a  solution  of  Eq.  (1(5)  to  within  about 
2  percent  for  all  altitude’s  above  104  km. 


ENERGY  (KeV) 


Figure  6.  Incident  Maxwellian  Energy  Distributions 
for  Characteristic  Energies  of  5,  2,  and  1  KeV  Each 
Containing  1  erg/cm2s 


Figure  7.  Energy  Deposition  Rates  from  the 
Analytic,  Range,  and  Strickland  Models  for 
the  5,  2,  and  1  KeV  Isotropic -Maxwellian 
Sources  as  in  Figure  (> 


We  are  interested  in  knowing  how  well  the  analytic  model  does  on  the  electron 
flux  itself.  Since  the  model  contains  the  forward  scattering  approximation,  we 
are  specifically  referring  to  *he  downward  moving  flux.  To  investigate  this  point, 
fluxes  from  the  analytic  and  Strickland  models  have  been  compared  at  selected 
altitudes  and  energies.  Figure  8  allows  for  such  a  comparison  at  110  km,  for 
the  2  KeV  isotropic-Maxwellian  case.  As  expected,  the  analytic  results  are  high 
near  m  ~  l  and  low  as  ju  ^  0  due  to  the  forward-scattering  approximation.  The 
differences  increase  with  decreasing  energy  since  scattering  becomes  more  im¬ 
portant  at  low  energy .  Figure  9  shows  the  hemispherically  averaged  flux  for  the 
two  models  from  ^  1  to  10  KeV  at  110  and  IGO  km. 


12.  CONCLUSION 

From  the  preceding  results  we  conclude  that  the  forward  scattering  and  aver¬ 
age  discrete  energy-loss  approximations,  which  at  first  seem  quite  severe,  to¬ 
gether  with  the  pseudoparticle  method  of  approximating  the  sums,  produce  a  sol¬ 
ution  for  the  primary  auroral  electron  flux  in  the  downward  hemisphere  which  is 
quite  simple  to  use  and  is  surprisingly  close  to  most  of  the  detailed  numerical 
results  of  Strickland  et  al.  The  formulae  we  give  for  the  primary  auroral 
electron  flux  and  the  quantities  derived  from  it  give  us  insight  into  the  physics  of 
the  electron  precipitation  process  and  also  provide  a  means  of  making  rapid  cal¬ 
culations  for  the  purpose  of  analyzing  auroral  data.  We  also  hope  that  the  reader 
will  find  he  results  presented  in  this  paper  useful  in  further  auroral  studies. 
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Figure  8.  Differential  Fleet  ron  Fluxes  from  the  Analytic  (N  •  12)  and  Strickland 
Models  at  110  km  for  the  2  KeV  Isotropic -Maxwellian  Source 


2  KeV  MAXWELLIAN 


Figure  9.  Ilemi spherically  Averaged  Fluxes 
from  the  Analytic  (N  =  12)  and  Strickland 
Models  at  110  and  160  km  for  the  2  KeV  Iso¬ 
tropic -Maxwellian  Source 
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